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C^) , Abstract 

In this work we present analytic and numerical treatments of the grav- 
itational wave signal from a pulsar which includes spindown. We consider 
Oh' phase corrections to a received monochromatic signal due to rotational 

and elliptical orbital motion of the Earth, as well as perturbations due to 
Jupiter and the Moon. We discuss the Fourier transform of such a signal, 
' which is expressed in terms of well known special functions and lends itself 

to a tractable numerical analysis. 
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(N ! 1 Introduction 

\ The detection of gravitational waves (GW) from astrophysical sources is one 

n ■ of the outstanding problems in experimental gravitation today. Gravitational 

ly-j ; wave detectors like the LIGO, VIRGO, LISA, TAMA 300, GEO 600 and AIGO 

are opening a new window for the study of a great variety of nonlinear curvature 
phenomena. Detection of GW necessitates sufficiently long observation periods 
O to attain adequate Signal to Noise ratio. The data analysis for continuous GW, 

' for example from rapidly spinning neutron stars, is an important problem for 

ground based detectors that demands analytic, computational and experimental 
to '. ingenuity. 

■ In recent works ^ we have implemented the Fourier transform (FT) of 

the Doppler shifted GW signal from a pulsar with the Plane Wave Expansion in 
Spherical Harmonics (PWESH). It turns out that the consequent analysis of the 
rS | Fourier Transform (FT) of the GW signal from a pulsar has a very interesting 

and convenient development in terms of the resulting spherical Bessel, general- 
ized hypergeometric, Gamma and Legendre functions. These works considered 
frequency modulation of a GW signal due to rotational and circular orbital mo- 
tions of the detector on the Earth. In the present analysis, rotational and orbital 
eccentric motions of the Earth, as well as perturbations due to Jupiter and the 
Moon, and pulsar spindown are considered. This formalism has a nice analytic 
representation for the GW signal in terms of the special functions above. The 
signal can then be studied as a function of various parameters associated with 
the GW pulsar signal, as well as the orbital and rotational parameters of the 
Earth. A brief analysis considering the parameters for the Lunar and Jovian 
perturbations is included. A detailed analysis of the spindown and perturba- 
tion effects has been recently done 3 . The numerical analysis of this analytical 
expression for the signal offers a challenge for efficient and fast numerical and 
parallel computation. 
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In this work, we present a variation of a parametrized model of pulsar spin- 
down discussed by Brady et al. |Jj and Jaranowski et al. [5] for the GW 
frequency and the phase measured at the ground based detector that includes 
the spindown parameters. The various types of motion of the detector can 
be naturally incorporated into this parametrized spindown model. Effects of 
amplitude modulation (AM), though not considered in this paper, can also be 
incorporated into this formalism without great difficulty We discuss the Fourier 
transform of the GW signal including the orbital and rotational corrections. In 
addition, a numerical study of the effects of the various detector motions on the 
phase of the GW signal and the dependence of the received signal on detector 
position and pulsar direction have been done. This might be relevant to GW 
detectors. 



2 Gravitational Wave Signal with Spindown Cor- 
rections 

Starting from a slightly modified form of the parameterization given by Brady 
and Creighton (2000)0), we obtain the following expression for the phase of a 
gravitational wave signal [H] : 



(t> A) = tt/o 
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The T k {z) (-1 < T fe (z) < 1, k > 0, -1 < z < 1) are the Chebyshev 
polynomials, and x = t + - • h, where f(t) is the position vector of the detector 
in the Solar System Barycentre (SSB) frame, and n is a unit vector in the 
direction of the pulsar. r m i„ is the spindown age of the pulsar in years, and /o 
is the GW frequency without spindown. The series in this expression sums to 




We have plotted this expression as a function of x for T. min = 40 and 1000 
years for z — —0.8, 0.7 and 0.9 (see Appendix, Figure 1). 
The GW signal is now given by 
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We have used the binomial expansion inside the square brackets. We define 
the Spindown Moment Integrals as follows: 



e-M*+^) . e - i2 -/' . ( X+ 2 ) A, (4) 
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then the fcth Spindown Moment Integral can be written as 

1 d k I n 
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Thus the Fourier transform of the signal can be written in terms of the partial 
derivatives with respect to fo of I generic-, which can be evaluated analytically. 

3 Corrections to the Position Vector of the 
Detector 

In this section we briefly outline the corrections to the position vector r(t) to 
account for the Keplerian ellipse, Earth's rotation, and the perturbations due 
to Jupiter and the Moon. 

An expansion of the elliptical orbit of the Earth, to second-order in eccen- 
tricity, leads to the following expressions for the x and y components of R or b, 
the vector specifying the Earth's orbital position (the z component is 0) : 



x(t) = a I 1 — ^e 2 ) cosAf + — cos2Af + ^ae 2 cos3Af — ^ae (7-a) 
8 / 2 8 2 



y(t) = a (l - ^e 2 j sin M + y sin2Af + ^ae 2 sin3AI - ^ae 2 sin M (7-b) 

(a = Sun-Earth distance, M — aj or f,i, e = eccentricity of orbit.) If the radius of 
Jupiter's orbit is much greater than that of the Earth, we can consider the Sun, 
to a first approximation, as moving around the Sun- Jupiter (S-J) barycentre, 
and the Earth orbiting the Sun. Then, the perturbation Rj due to Jupiter is 

m 

Rj = [Rj (cos (wjt) - 1) , Rj sin (tujt) , 0] (8) 

where u>j is the orbital angular frequency of Jupiter, and Rj is the distance of 
the Sun from S-J barycentre. 
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Similarly, the vector from the Earth-Moon barycentre to the Earth is: 



Rm 



Rem (cos - 1) 

Rem sin (w M t) 




0) 



Here u>m is the Moon's sidereal orbital angular frequency, and Rem is the 
distance from Earth to the Earth-Moon barycentre. Also, the vector specifying 
the detector's position due to the rotating Earth, R ro t, is given by 



Rrot — 



Re sin a (cos (u> ro tt) — 1) 
Re sin a sin (uj rot t) cos e 
Re sin a sin ((aj rot t) sine 



(10) 



where Re = Earth's radius, a — detector co-latitude, e = angular tilt of Earth's 
axis, and Lo ro t = Earth's sidereal rotational angular frequency. The position of 
the detector is now given in the form: 



r(t) = R orh {t) + R rot {t) + Rj(t) + R M (t) 



(11) 



4 Contributions of Perturbations to the Phase 
of the GW Signal 

The contribution of a perturbation to the phase of the GW signal (in radians) 
is given by </>(t) = 2Trfo r '^ t < }' n , where fi(t) is the contributing vector for the per- 
turbation |2 • We have plotted the phase contributions for the Earth's circular 
and elliptical orbital and rotational perturbations, as well as contributions from 
the Jovian and Lunar perturbations, as functions of 9 (the angle of the pulsar 
from the orbital plane normal) and <fi (the angle of the pulsar from the £-axis in 
the orbital plane), for t = 6 lunar months (6 x 29.5 d), a = 5 and /o = 1 kHz 
(see Appendix, Figure 2.) In each plot, the maximum phase shift occurs near 
6 = | (when the pulsar direction is in Earth's orbital plane). The Rvalue at 
the maximum is more variable. The approximate maximum phase shifts for the 
above perturbations are, respectively, 6.3 x 10 6 , 3.2 x 10 6 , 1.8 x 10 2 , 2.6 x 10 4 
and 1.2 x 10 2 radians. For comparison, the phase of the signal in this time 
period without Doppler shifts (27r/of) is 9.6 x 10 10 radians. 

We have also plotted the positional errors from using the circular approxi- 
mation versus those of equation 7, for a 1 kHz gravitational wave (see Appendix, 
Figure 3). Four orders of magnitude are gained, a substantial reduction in the 
error. In addition, we see that the phase error can be quite large if we ignore 
the effects of Jupiter and the Moon over long integration times. 
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5 The FT of the Perturbation-Corrected GW 
Pulsar Signal 

In Valluri et al. (2005) [3], we arrive at an analytic expression for the Fourier 
transform of a GW signal, considering the perturbations above: 
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Here c is the speed of light, and ljq = 2tt/q. Also, 
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Feynman has thoroughly discussed the diffraction pattern due to the factor 
in the curly brackets in the expression for ?/>2 which gives the resultant amplitude 
due to R equal oscillators (R — 365 in our case) . This factor is analogous to 
the expression for the scattering amplitude of an electron in a crystal lattice. 
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The intensity of the signal, 



Kf) 



can be expressed in analytic form from 
Equation 12. A detailed numerical analysis of this expression for the variety of 
parameters, which can take huge values (> 50000), present in the GW signal is 
a challenge for high performance parallel computation 8 . 



6 Conclusions 

We have presented in this paper the rudiments of a simple analysis for spindown 
for sources of continuous GW. For the more computationally intensive search 
over all sky positions and spindown parameters, it is important to be able to 
calculate the smallest number of independent parameter values which must be 
sampled in order to cover the entire space of signals. The PWESH improves 
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the numerical accuracy and convergence of analytic FT's, spindown corrections 
associated with the GW signal, and also enables the estimation of parameters 
and distributions relevant to GW Data Analysis. It has also found many im- 
portant applications, for example in the study of the multipole moments in the 
Cosmic Microwave Background Anisotropies [§]. The proper blend of analytic 
and numerical integration for accurate GW data analysis remains an interesting 
question to be explored further. 
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8 Appendix 
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Figure 1: Plot of summation as a function of x for r m « n = 40 years (left) and 
r mm = 1000 years (right) for different z values. 




(d) (e) 

Figure 2: Plot of phase shift due to (a) circular orbit, (b) elliptical orbit, (c) 
Earth's rotation, (d) Jupiter, and (e) the Moon, as a function of 9 and <p. 
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Figure 3: The difference between the Earth's position as computed by a circular 
approximation versus a Keplerian one, measured in cycles of a 1 kHz gravita- 
tional wave over the course of one year (dashed grey line) . The (much reduced) 
difference when equation 7 is used is shown as a solid black line. Approximate 
contributions to the motion of the Earth due to the Moon and Jupiter are shown 
for comparison. 



8 



